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ABSTRACT 



Dwarf galaxies provide insights on the processes of star formation and chemical enrichment at the low end of the galaxy mass 
function, as well as on the clustering of dark matter on small scales. In studies of Local Group dwarf galaxies, spectroscopic samples 
of individual stars are used to derive the internal kinematics and abundance properties of these galaxies. It is therefore important to 
clean these samples from Milky Way stars, not related to the dwarf galaxy, since they can contaminate the analysis of the properties 
of these objects. Here we introduce a new diagnostic for separating Milky Way contaminant stars - that mainly constitute of dwarf 
stars - and red giant branch stars targeted in dwarf galaxies. As discriminator we use the trends in the equivalent width of the nIR 
Mglline at 8806.8 A as a function of the equivalent width of Call triplet lines. This method is particularly useful for works dealing 
with multi-object intermediate resolution spectroscopy focusing in the region of the nIR Call triplet. We use synthetic spectra to 
explore how the equivalent width of these lines changes for stars with different properties (gravity, effective temperature, metallicity) 
and find that a discrimination among giants above the horizontal branch and dwarfs can be made with this method at [Fe/H]> -2 dex. 
For -2 < [Fe/H] < -1, this method is also valid to discriminate dwarfs and giants down to approximately one magnitude below the 
horizontal branch. Using a foreground model we make predictions on the use of this new discrimination method for nearby dwarf 
spheroidal galaxies, including the ultra-faints. We subsequently use VLT/FLAMES data for the Sextans, Sculptor and Fornax dSphs 
to verify the predicted theoretical trends. 
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1. Introduction 

There is much to be learned from the galaxies within the Local 
Group which can be studied in better detail than any other sys- 
tem, on a star-by-star basis. Within the Milky Way (MW) halo 
itself, we find satellite galaxies which span a range over 10 mag- 
nitudes in brightness. The great majority of the MW satellites 
belong to the class of early type dwarf galaxies, which are the 
smallest and least luminous galaxies known to date. These in- 
clude the ultra faint dwarfs (UFDs) which were relatively re- 
cently discovered with the Sloan Digital S ky Survey (SPSS) . 
Several future s urveys, like Pan- STARRS dKaiser et al.l 12002) 
and Skymapper dKeller et al.ll2007l) . are expected to further in- 
crease the number of detected dwarf galaxies around the MW, 
towards lower surface brightnesses. 

Great effo rt is devoted to the study of these small early type 
galaxies (see iTolstov et al.ll2009l for a review) as they provide 
precious insights on the processes that drive galaxy evolution at 
the low end of the galaxy mass function, as well as on the clus - 
tering of dark matter on small scales (e.g. lGilmore et al. 2007). 
Since stars in these galaxies are resolved, the internal proper- 
ties for the galaxies are derived from photometric and spectro- 
scopic studies of individual objects. This requires the application 
of methods to discard from the sample foreground/background 



* Based on FLAMES observations collected at the ESO, proposals 
171.B-0588, 076.B-0391, 079.B-0435 



stars which are not bound to the dwarf galaxy, but that are in- 
stead within the MW itself (hereafter "contaminants" or "inter- 
lopers"), and that would influence the properties derived for the 
dwarf galaxy. 

General methods to exclude obvious interlopers include the 
use of color-magnitude diagrams (CMDs) to select a region in 
magnitude and color where one would expect member stars in 
a certain stellar evolutionary stage to lie given the distance to 
the object. For most of the studies of Local Group early type 
dwarf galaxies a photometric selection box is placed on the red 
giant branch (RGB) - the most luminous feature in the CMD 
of ancient/intermediate age stellar populations. However, such a 
selection box will also contain contaminants, these consist pre- 
dominantly of much intrinsically fainter but close-by main se- 
quence dwarf stars from the MW thin or thick disk. In order 
to distinguish the contaminant dwarf stars from the RGB stars 
members to the dwarf galaxy additional discriminators need to 
be used. 

Spectroscopic studies can rely on a further selection cri- 
terion to weed out foreground contamination, by using the 
line-of-sight velocity of the individual stars as compared to 
the systemic velocity of the dwarf galaxy. This is commonly 
done in a statistical sense, either using sigma-clipping proce- 
dures (which imply a negligible probability to find a member 
beyond a certain velocity threshold) or maximum likelihood 
analyses, with more detailed implementations modeling both 
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the expected population of dwarf galaxy stars and MW stars 
(e.g. [Battaglia et alJl2008at IWalker et aill2009bt iMartinez et al l 
1201 ll) . Additional properties than line-of-sight velocities can be 
included in the analysis, e.g. distan ce of the star from the galaxy 
center (e.g. Battag lia et al.ll20 08a). and information on spectral 
indices dWalker et al.l2009bl for the Mg index around A 5 170 A ) 
or metallicity dMartinez et al.ll201ll) of the star. The most sophis- 
ticated of th ese implementations from a statistical po i nt of view 
are those of I Walker etail (l2009bl) and lMartinez et alJ (1201 ll) . 

However, even the most efficient of these methods cannot by 
construction select out MW contaminants with similar proper- 
ties (colors, apparent magnitudes, line-of-sight velocities etc) as 
the member stars to the observed galaxy. Depending on the type 
of stars surveyed, the systemic velocity of the galaxy and its po- 
sition on the sky (most importantly Galactic latitude), contami- 
nant stars satisfying all these requirements might still be numer- 
ous or at least statistically significant in number. Furthermore, a 
selection criterion based on position, velocity (and/or metallic- 
ity) needs assumptions about the underlying distribution and can 
therefore also bias the final analysis. For example, it is unclear 
whether features in the outer parts could be mistaken for con- 
taminants: if looking for tidally disrupted stars, one expects to 
find them in the very outer parts of the dwarf and at velo cities 
rather discrepant from the systemic (e.g. lMunoz et al.ll2006l for 
an observational study of tidal disruption in Carina). In statis- 
tical procedures those stars would most likely be considered as 
non- members. Similarly, sta rs in the outer "shell" of the Fornax 
dSph dColeman et al.ll2005h . arguably the remnant of an accre- 
tion event, would most likely to rejected because well outside 
the limiting radius of the dwarf as inferred by its surface bright- 
ness profile. Also, stars with very low and high metallicity with 
respect to the rest of the population, could be rejected as non- 
members, while true members with deviating metallicities could 
give important insights on the earliest phases of chemical en- 
richment and on the capability of the dwarf to retain gas and/or 
metals, respectively. 

These examples illustrate the benefits of having also methods 
that weed out contaminants on the basis of sensitivity to gravity, 
since the stars targeted in external galaxies are RGB stars and the 
foreground consists mostly of dwarf stars from the thin or thick 
disk. Such a discriminator can be used separately or combined 
with other (statistical) methods to increase the efficiency of the 
elimination of interlopers. 

Photometric techniques are present in the literature which al- 
low a distinction between dwarf and giant stars by constructing 
indices sensitive to gravity and luminosity from combinations of 
particular filters. For example, the DD051 filter is centered on 
the Mg b/MgH feature near 5 170 A , which is strong in late-type 
dwarfs and weak in giants (e.g. Morrison et al. 2001). Several 
works in the literature have shown the effectiveness of separat- 
ing giant an d dwarf stars u sing a c ombi nation of Washington and 
DDO filters: iMunoz et aD (see e.g. l2006l) pre-selected photomet- 
rically their spectroscopic targets using the M51 filter, finding 
probable members to the Carina dSph out to very large distances 
from its center, where genui ne members would have otherwise 
been buried in the foreground. Fari a et all d2007l) used Stromgren 
photometry on the Draco dSph, showing that the Stromgren c\ 
index, sensitive to luminosity, can be used together to the (b — v)o 
color to separate stars in various evolutionary stages. While pho- 
tometric techniques provide broad information on the star for- 
mation history of the dwarf galaxy a nd in some cases on the 
metallicity distribution of its stars (e.g. Faria et al. 2007), infor- 
mation on the line-of-sight velocity of the individual stars ap- 



pears crucial to assess membershi p in some case s even when 
using Washington photometry (e.g. lMorrison et al. 200lh. 

There also exist spectroscopic studies of early type Local 
Group dwarf galaxies that have used features in the spec- 
tra of the stars to dis crim inate between giant and dwarfs. 
Spinrad & Taylor! d!971l) and lSchiavon et all d 19971) showed the 
Na I doublet of absorption lines at 8183, 8195 A to be strongly 
dependent on gravity and temperature. Several studies have 
adopted the equivalent width (EW) of the Na I doublet as a 
diagnostic to discriminate between dwarf and giant stars, with 
the dwarf stars displaying much larger EWs than giant stars ( see 
iGilbert et al.1l200a ISimon & Gehall2007l: iMartin et alJl2007l on 
the halo of M3 1 for t he former, and on U FDs for the latter two 
works). In their work. lGiibert et al] d2006l) showed that this diag- 
nostic offers clear differentiation at colors (V-I)> 2, while there 
appears to be significant overlap in the distribution of Na I dou- 
blet EW of dwarf and giant stars for bluer colors. 

The adoption of one, or more, method to assess membership 
depends on the particular scientific problem, galaxy, and data- 
set under consideration, and the fraction of contaminants that 
can be tolerated so as not to significantly affect the results. In 
some cases the presence of interlopers has a dramatic impact on 
the conclusions. For example mass estimates derived from mea- 
surements of the line-of-sight velocity dispersion of samples of 
individual stars in LG early type dwarf galaxies have been used 
to discuss the poss i ble existence of "un iversal" mass value (e.g. 
Strig ari et al.l2008tlWalker et al.l2009al) . In their study of UFDs, 
ISimon & Gehal d2007l) adopted a wealth of membership indica- 
tors, such as spatial position of the star, location on the CMD, 
EW of the Na I doublet, line-of-sight velocity and metallicity, 
and found for the Hercules dSph an internal velocity dispersion 
of 5.1+ 0.9 km/s. By remov ing MW stars using Stromgren pho- 
tometrv. lAden et all d2009al) found a lower velocity dispersion in 
their cleaned sample , 3.7+0.9 km s" 1 . In this particular example, 
the derived mass of lAden et al.l d2009bl) from this new velocity 
dispersion places Hercules well bel ow the "common m ass scale" 
for MW early type dwarf galaxies dStrigari et al.ll2008l and ref- 
erences therein). 

Another instructive example is the case of Willman 1, whose 
apparent spread in the metal abundance distribution was used in 
favor of the classific ation as a very low luminosity dwarf galaxy 
(Mart in et al.l 12007). Subsequent studies based on high resolu- 
tion spectroscopic observations revealed that a large percentage 
of stars classified as members on the basis of thei r position on 
the CM D, velocity and EW of the Na I doublet bv lMartin et al.l 
(120071) were either foreground dwarfs or halo giants, renewing 
the possibility that Willman 1 could be a disrupted metal-poor 
globular cluster. 

The examples above show the importance of careful inter- 
loper removal by using as many membership indicators as possi- 
ble. This is particularly important as the galaxies studied become 
fainter and consist of less stars, and for studies which explore 
the outer regions of dwarf galaxies, where the expected ratio of 
members/interlopers is large. 

In this paper we present an additional discriminator to sep- 
arate dwarf and giant stars, i.e. the EW of the Mg I line at 
8806.8 A used in combination with the sum of the EW of the 
two strongest Call triplet (CaT) fines (£W), at A = 8542.1 and 
8662.1 A (see Figure [p. This method was empirically intro- 
duced by Battaglia et alJ (1201 ll hereafter Bll) as a criterion to 
discriminate RGB stars belonging to the Sextans dSph from the 
MW contaminant stars. Using a data-set with a limited extent 
in metallicity, the authors noted that the great majority of the 
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stars with line-of-sight velocities very close to the systemic of 
Sextans (i.e. highly likely member RGB stars) clustered around 
small values of the Mg I EW, while stars with velocities more 
than 4<x away from the systemic (i.e. highly likely non-members) 
clustered at larger values, > 0.5 A . In this paper we extend the 
empirical work of B 1 1 : we use synthetic spectra to explore the 
validity of the method over a range of metallicities, gravities and 
£W, demonstrating the effectiveness of using the Mg I line in 
combination to the CaT lines in removing MW contaminants. 

This method is a convenient choice: Mgl and CaT lines 
are close in wavelength space and relatively broad, which 
means they will typically both be present and measurable in 
intermediate-resolution studies of this wavelength regime. They 
are located in the red part of the spectrum, where the target 
RGB stars are brightest, with optimization of the exposure time. 
Furthermore, the wavelength region including CaT lines has be- 
come increasingly more used in intermediate resolution spectro- 
scopic studies of LG early type dwarf galaxies: CaT lines are 
both very suitable for radial velocity determinations and offer 
the possibility of deriving estimates of [ Fe/H] accurate to ±0-15- 
0.2 dex over a wide metallicity range (Starkenb urg et alJ feOlO) 
for large numbers of individual stars and much lower observing 
times with respect to high resolution studies (typically a S/N~ 10 
is required). The method can be applied to a wealth of galaxies 
observed in this wavelength regime. In the ESO archive a large 
body of data already exists at intermediate resolution in the only 
set-up including the CaT, which covers the region between 8206 
and 9400 A. This is by far the most frequently used set-up for the 
FLAMES instrument. Note that the Na I doublet is just outside 
this wavelength range. The current archive includes observations 
of Fornax, Sculptor, Sextans, Carina, Leo II, Leo IV, Hercules, 
Boo I, IC 1613, Sag, Segue 1 in this wavelength range, by var- 
ious groups. The CaT region has also been targeted by stud- 
ies using other telescopes than VLT, such as the work of Koch 
et al. 2007 on Leo I. Furthermore, the future large GAIA-ESO 
Survey (GES) will target the wavelength region covering the nIR 
CaT and Mg I line (V. Hill, GES Team, private communication). 
Several existing and future data-set could therefore benefit from 
this discriminator. 

The paper is structured as follows. In Sect. [2] we use syn- 
thetic spectra to explore the trends between the E W of Mg I line 
and CaT lines for stars of different metallicity and gravity and 
to derive a relation that allows us to discriminate among giant 
and dwarf stars on the basis of these lines. We then comment 
on the expected success rate of this method along the line-of- 
sight to different MW satellites in Sect. [3] In Sect. |4] we use 
data for the Sextans, Sculptor and Fornax dSphs fr om the Dwarf 
Abun dances and Radial Velocities Team (DART, Tolsto v~et al.l 
2006), to explore the behavior of the lines empirically for galax- 
ies with different characteristics, compare the results to the mod- 
els and show the agreement between models and observations. 
Finally we apply the method to the DART data-set in Sect. 14.11 
and present our conclusions in Sect. [5] 



2. The lines and their sensitivity 

In order to investigate the behavior of the nIR CaT and 
Mg I absorption lines we first have to define the boundaries of 
the physical parameter space we are working in. In practise, 
we do this using isochrones covering possible ranges of [Fe/H] 
and ages as covered by the stars in Local Group dwarf galax- 
ies and MW foreground in order to select the appropriate re- 
gions in log(g) - T e ff space (see Figure|2]i. The isochrones we use 
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Fig. 1. VLT/FLAMES spectrum from DART observations of a 
star along the line-of-sight to the Sculptor dSph (heliocentric ve- 
locity = 122.7 km s _1 and signal-to-noise ratio S/N/A =64). The 
lines used in this work, i.e. the two strongest CaT lines and one 
Mg I line, are marked. 



are Yonsei-Yale isochrones (e.g. lYi et al.l l200H lDemarciue et al.l 
2004) for 1, 5 and 12 Gyr, covering the wide range -3.6 < 
[Fe/H] < +0.0. To represent general trends within the MW stars, 
the alpha enhancements chosen are +0.4 for [Fe/H]< -1 dex 
and +0.0 for the solar [Fe/H]. However, since spectroscopi c 
observatio ns of stars in th e Sculptor ( e.g. iTolstov et al. l 12009). 
Fornax dLetarte et al.l2010h and Carina dKoch et al.l2008l) dSphs 

show [cc/Fe] around solar values already at [Fe/H] 1 dex, for 

completeness we add an isochrone set of [a/Fe]=0.0 dex and 
[Fe/H]= -1 dex. 

Because the dwarf galaxies currently studied in large spec- 
troscopic surveys ( e.g I Walker et al.ll2006h ISimon & Gehall 2007: 
lAden et al.l l2009al: IWalker et alJ l2009rA iKirbv et alJ 1201 Oh are 
distant objects, typically only the giants above the horizon- 
tal branch (HB) will be bright enough to be followed up with 
spectroscopy; in systems with sparsely populated RGBs though 
observations often need to be p ushed to lower magnitudes, 
theref ore below the HB (e.g. lMartin et al.ll2007l: ISimon & Gehal 
2007). The isochrones plotted in Fig. [2] show that giant stars of 
any metallicity or a ge above the HB (which will have <My < 
+ 1 tvpicallv lRich et al.ll2005trCatelan & Cort es 2008), and most 
of those below the HB, will not be found at bluer colours 
than V-I ~0.85, correspond ing to stars hotter than ~5250 K 
dRamfrez & Melend ez 2005). In our following analysis, aiming 
at disentangling the populations of giant and dwarf stars which 
would be pre-selected photometrically, we therefore will not take 
temperatures > 5250 K into account. This temperature selection 
cuts out an important fraction of the hot, halo and thick disk 
stars in the foreground contamination. Figure [2] also shows that 
at the lowest metallicity end the isochrones become less sensi- 
tive to metallicity and thus we expect no significant changes if 
even lower metallicities ([Fe/H] =-4 dex and below) would be 
considered in the analysis. 

On the other end of the temperature scale we place a cut-off 
at V-I= 2.1, corresponding to a temperature of 3600K. Although 
in theory very metal-rich giants stars can reach even cooler tem- 
peratures, these are rare in dwarf galaxies because most dwarf 
galaxies have low metallicity. Indeed, following this consid- 
eration, most photometric boxes that have been used for pre- 
selection of spectroscopic targets are well within this range. Note 
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that with this cut, the only population of giant stars completely 
excluded from the analysis is the one of super-solar metallicity 
and with an age of 12 Gyr, unexpected in such systems. 

Since we address in this work the validity of the use of 
the Mgland CaT lines to discriminate RGB stars in Local 
Group dwarf galaxies from MW contaminants, we use in the 
analysis synthetic spectra with physical properties (log(g), T e ff, 
etc) covering the range expected for both populations of stars 
within the color range discussed above. We have calculated a 
grid of synthetic spectra for a wide range of metallicities, tem- 
perat ures and gravity using (OS)MARCS model atmospheres 
(e.g. Gustafsson et al. 2008; Plez 2008) and the Turbospectrum 
program dAlvarez & Plezl[l998h . updated consistently with the 
iGustafsson et alJ d2008l) MARCS release. In agreem ent with the 
trend s followed by most of the stars in the MW (e.g. lVenn et al.l 
120041) . we adopt MARCS models with "standard" chemical 
composition, i.e. the values chosen for the alpha elements (taken 
as O, Ne, Mg, Si, S, Ar, Ca, and Ti) relative to Fe for the mod- 
els are [a/Fe]=+0.4 at the metallicities below solar, linearly de- 
creasing towards [a/Fe]=0.0 at solar [Fe/H]. The microturbulent 
velocities used in our mo dels are 1-2 km s" 1 , in agreement with 
the observed range (e.g. Bar klem et al.ll2005l) . 

To allow a comparison of the data with the synthetic spec- 
tra we define regions in T e ff-log(g) space where the giants 
above the HB, (sub)giants below the HB and main-sequence 
branches would lie (see upper panels of Figure [3]). For "clas- 
sical" dSphs, which have a well populated RGB, spectroscopic 
samples mostly target RGB stars above the HB, therefore the 
member stars we are looking for will all fall within the 'giant' 
boxes (green boxes in Figure [3). For those UFDs for which ob- 
servations must be pushed below the HB, members will also fall 
within the 'sub-giant' boxes. For all the models that fall within 
either one of the selection boxes we measure the EW of the 
Mg I line and of the two strongest CaT lines in the correspond- 
ing synthetic spectrum. We measure the EW of the CaT lines 
using a Gaussian fit to the individual lines, corrected by a fac- 
tor 1.1 to include the non-Gaussian wings. This factor is deter- 
mined by comparing the EW from the Gaussian fit to the one 
from the integrated line measure (f or mor e details see Battaglia 
et al. I2008bl and Starkenburg et al. l2010h . We also correct the 
EWs of the CaT lines for non-LT E effects using the equations 
given in Starkenbu rg et al.l ((2010), these corrections are mainly 
important for the very low-metallicity regime. The Mglline is 
measured using a simple integration of the area under the con- 
tinuum of 6 A around the line, as also done for the data (see 
Sect.gji. 

In the bottom panel of Figure[3] we show the sum of the EW 
for the two strongest CaT lines (SWcat - ) versus the EW of the 
Mglline (EWM g ). For [Fe/H] = -2 and -1 dex the dwarf stars 
can be separated from the giants and sub-giants. At solar metal- 
licity, the distinction remains clear between the dwarfs and the 
giants, but not for the sub-giants. Note though that, given the 
metallicity-luminosity relation (see iKirbv et al]l2010l for early 
type dwarf galaxies in the Local Group), stars with metallicities 
close to solar are not expected in UFDs, which are the systems 
for which one often needs to target stars below the HB due to 
their sparsely populated RGB. The figure also shows that the 
distinction among the different classes of objects would not be 
possible if one were to collapse all data onto one axis (i.e. would 
use just one of the lines): it is the combination of these lines 
that makes the distinction between the giants and dwarfs possi- 
ble in most cases. In the extremely low-metallicity regime (i.e., 
[Fe/H] = -3), giants, sub-giants and dwarfs appear to occupy the 
same region in the EWMg, £Wc 2 t plane. However this is not a 




Fig. 2. Colour- magnitude diagram using Yonsei- Yale isochrones 
for 1 Gyr (dashed lines), 5 Gyr (dotted lines), and 12 Gyr 
(full lines) and metallicities [Fe/H]=-3.6 (lowest value possi- 
ble within the Yonsei-Yale set in black), [Fe/H]=-3.0 (grey), 
[Fe/H]=-2.0 (red), [Fe/H]=-1.0 (blue), [Fe/H]=0.0 (green) and 
[Fe/H] =+0.5 (purple). [a/Fe] is set to +0.4 for all isochrones 
except for [Fe/H]=0.0 and [Fe/H]=+0.5 which have [a/Fe]=0.0 
and an extra set of isochrones plotted for [Fe/H]=-1 and 
[a/Fe]=0.0 (cyan). The vertical line roughly indicates the onset 
of the RGB. 

concern because in Sect. [3]- where we discuss the expected con- 
tamination for different objects - we will show that interloper 
stars with such low metallicities are very rare. 

A distinction line between giants and dwarfs is overplotted 
in all panels and is: 



EW Mg (mA) = 



300 



ifA< 3750 mA 



0.26 x EW CaT - 670.6 if A > 3750 mA 



(1) 



As we will discuss in Sect. [3] the fraction and type of con- 
taminant stars, and therefore their location and importance on the 
EWm s vs SWcaT plane, will change according to the specifics 
of the dwarf galaxy (position on the sky, systemic velocity, dis- 
tance). This line therefore provides a criterion to keep in as many 
of the member stars of the dwarf galaxy, and throw away as many 
contaminants as one can. 



3. The foreground contamination 

3.1. Fornax, Sculptor and Sextans 

As a test case of different mixtures of MW interlopers we look at 
the expected contamination along three different line-of-sights, 
specifically to the Sextans, Sculptor and Fornax dSphs. 

To study the likely foreground cont amination towards these 
objects, we used the Besancon model (Robin et al. 2003]). All 
interloper stars were selected within a 4 deg 2 solid angle around 
these galaxies, at a distance between and 200 kpc. As in Sect.|2] 
we restrict the analysis to the color range 0.85 <V-I< 2.1. We 
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Fig. 3. Upper panels: Yonsei-Yale isochrones for 1 Gyr (dashed lines), 5 Gyr (dotted lines), and 12 Gyr (full lines) for one particular 
metallicity. For [Fe/H]=-3 and [Fe/H]=-2, [a/Fe]=04, for [Fe/H]=0.0, [a/Fe]=0.0 and for [Fe/H] both alpha enhancements are 
given ([a/Fe]=0.0 in grey). All values for which the synthetic spectra are calculated are covered with a grey dot. Overplotted are 
regions where the giants above the HB (green), (sub)giants below the HB (blue) and main-sequence (red) branches would lie in 
Teff-log(g) space, if the temperature restrictions from photometry are taken into account (vertical dotted lines). Lower panels: EWs 
for the two strongest CaT line versus the Mgl 8806.8 line for models that fall within the boxes of the panel above. Giants above the 
HB are shown as green triangles, (sub)giants below the HB as blue diamonds and dwarfs as red asterisks. Overplotted a distinction 
line for dwarfs and giants (except for the extremely low-metallicity case), this line is the same in all panels. 



also make a selection in apparent magnitude, which would in- 
clude giant stars above the HB at the distance of these "classical" 
MW dSphs. For this we use the distance moduli from the liter- 
ature and information on the HB magnitudes (see Table [U, and 
we adopt an absolute magnitude in /-band M/ = -4 for the tip 
of the RGB. This range in apparent magnitude will be populated 
also by a range of MW contaminants - mainly main sequence 
dwarf stars - with the appropriate combination of absolute mag- 
nitude and distance from the Sun. 

The leftmost panels of Figure [4] show the foreground sam- 
ple in the direction of these three dwarf galaxies (note that in all 
these panels the contour levels are not linear), while the number 
of contaminants is given in Table Q] It is clear that the type of 
stars within the foreground sample will depend on the direction 
in the sky. Sculptor, being the galaxy at highest latitude has the 
least foreground in general, but a relatively high fraction of its 
foreground are giant stars, difficult to weed out with our method. 
Although all these galaxies are at high(ish) latitudes, there still is 
a very significant contamination from the disk in all cases, in par- 
ticular from the thick disk. The thin and thick disk contribute al- 
most entirely dwarf stars (log(g) > 4.3), while the fraction drops 
to about 50% for the contribution of dwarf stars from the halo. 
As previously discussed, the figure also shows that the great ma- 
jority (>90%) of the contaminants are found at [Fe/H]> -2 dex, 
where the distinction between giants and dwarfs is easier. The 



stars with [Fe/H]< -2 dex mostly belong to the stellar halo, 
which is the least dominant population in all of these cases. 

In previous work, the member stars of the dwarf galaxy were 
selected out from the contaminants on the basis of their line-of- 
sight velocities. The stars were for example required to have a 
line-of-sight velocity within 3cr from the systemic velocity of 
the galaxy (for the Fornax dwarf galaxy 2.5cr was used, because 
of its systemic velocity is closer to the velocity of the Galactic 
disk[Battaglia et al. 2006). Using the simulated catalogs of MW 
contaminants we find that in all cases, as expected, some fore- 
ground stars will remain which have similar velocities as the 
dwarf galaxy stars, as is shown in the second column in Figure|4] 
The fraction though changes significantly from galaxy to galaxy. 
For Sculptor and Sextans the fraction of interlopers retained by 
the velocity criterion is low, as listed in TableQ] but the cleaning 
by velocity is much less effective in the Fornax dwarf galaxy, 
where 44% of the interloper stars have a velocity within 3-cr of 
the systemic velocity of Fornax (35% if 2.5-cr is used). In com- 
parison, for Sextans this is just 0.9% within 3-cr. 

In the third column of Figure |U we weed out contaminants 
according to the "line criterion", i.e. for which the closest syn- 
thetic spectrum in physical properties has a EWM^/SWcaT ra- 
tio which falls above the dividing line in Figure Q\ To derive 
the information on the line criterion all interloper stars from the 
Besancon file were linked to the closest model within the grid of 
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Table 1. Properties of the foreground contaminati on towards the MW satellites. The colum ns list: the name of the galaxy (1), its 
galactic coordinates (2,3); distance m odulus from | Mateol (1 19981) for t he classical satellites, iBelokurov et"al] (2007) for Coma B, 
Hercules, CVe nll, Leo IV and Segu e L lZuckeret al. ( 2006al) for UMall. lZucker et alj d2006bl) for CVenI an dlSimon & Gehal d2007h 
for UMal and IMartin et al.l (l2007h fo r Willman 1 (4); V-band magnitude of the horiz ontal branch from iMartin et all d2007l) for 
UMal and CVenI, Aden et al.1 d2009ah for Hercules and [Irwin & Ha tzidimitriou ( 1995) for the cla ssical satellit es (5); heliocentric 
systemic ve locities from B 1 1 for Sextans. IBattaglia et alj d2006l) for Fornax. IBattaglia et all d2008bl) for Sculptor. iKoch et al.l f2006) 
for Carin a. | Mateol (1 1998h f or the classical d warf galaxies and Simon & Gehal (120071) for most ultra-faint dwarf galaxies, except 
Segue I dSimop et al.ll2010h and W illman 1 dMartin et alj|2007l) (6); corresponding velocity dispers ion from the same references 
except we use Walker et al. (2009b) for the classical dwarf galaxies with systemic velocities from Mateo] d!998l) (7); the number 
of expected MW contaminants within a solid angle of 4 deg 2 along the line-of-sight to the galaxy (8); The last 4 columns list the 
percentage of contaminants retained in the sample when using only the velocity criterion (9), the line criterion (10), when applying 
both (11) and applying both relaxing the velocity selection to 4-cr (12). 
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synthetic models (described in Section |2]i in [Fe/H], log(g) and 
temperature space and the EWs of the CaT and Mg I lines of the 
closest model are used as their measurements. In all three galax- 
ies it can clearly be seen that although a lot of contaminants are 
selected out (only 12 - 17% remains), the foreground consisting 
of giants remains (see also Figure [5] for the distribution of con- 
taminants in the EWMg vs 2Wc 3 t plane before and after veloc- 
ity selection). The line criterion is less effective than the velocity 
criterion for Sextans and Sculptor. From Table[T]however is clear 
that for Fornax, that has a systemic velocity much closer to the 
one of MW thin disk with respect to the other two dSphs, the 
line criterion performs much better than the velocity criterion on 
its own. 

The best results in all cases are obtained when the velocity 
and the line criterion are used in addition to each other, as they 
are completely complementary in the sense that they will select 
out different sets of foreground stars. The line criterion will be 
selecting especially on the physical properties (mostly gravity) 
of the star, the velocity criterion just on its dynamics within the 
galaxy. Especially metal-rich dwarf foreground stars which have 
similar velocities and will be retained within the sample using a 
velocity cut, will be very easily separated by the method using 
absorption lines. In the last column of Figure [4] we show the 
result of a combination of both methods. In all galaxies a very 
large fraction of the interlopers can be weeded out and only a 
few percent of the interlopers still remain. 

An additional value of complementary use of both a velocity 
and line criterion is given in the last column of TableQ] where we 
relax the velocity criterion allowing for a selection of member 
stars within a larger range of velocities, e.g. 4-cr, and find that 
this does not cause a large increase in the number of interlopers. 



This can prove extremely useful, since one of the concerns of 
imposing strict velocity criteria is the possible biases introduced 
by the exclusion of member stars with more deviant velocities 
from the mean, for example when investigating the mass content 
of the galaxy from the velocity distribution of its stars. 

3.2. Predictions for the use of the method for other Milky 
Way dwarfs, including the ultra-faints 

We now extend the previous analysis to the other early type MW 
satellites. In the lower part of Table[T](below the horizontal line) 
we summarize the relevant information for all these galaxies for 
which dynamical analyses have been carried out, including the 
UFDs dStrigari et al.ll2008l) . For each of these galaxies we have 
extracted a catalog of foreground stars from the Besan^on mod- 
els and analyzed it in the exact same way as described for the 
classical satellites Sextans, Sculptor and Fornax. In the various 
galaxies spectroscopic samples extend down to different depths, 
and most UFDs have been surveyed in observations also below 
the HB to increase the number of targets due to their sparsely 
populated RGB. For uniformity with the previous analysis we 
consider here only the region above the HB; note though that this 
method will be effective at the metallicities expected for stars in 
UFDs also down to approximately one magnitude below the HB, 
corresponding to the sub-giant branch boxes in Fig. 3 (see also 
the Appendix for an exploration of the method performance to 
fainter magnitudes and bluer colors). 

As it can be seen from Table [1] in all galaxies the line crite- 
rion can weed out approximately 80-90% of the interlopers. The 
efficiency of a velocity criterion on the other hand, is much more 
dependent on the particular galaxy and its systemic velocity. In 
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Fig. 4. Left panels: [Fe/H] and log(g) for all foreground stars drawn from the Besancon model towards the direction of the dwarf 
spheroidal galaxies Sextans (top row), Sculptor (middle row) and Fornax (bottom row). Contours are shown on a 25 by 25 grid and 
represent percentages of 0.001, 0.01, 0. 1 and 0.5 the total number of stars. The rest of the panels show subsequently the distribution 
of foreground stars left after applying a velocity criterion of 3<x around the systemic velocity of the galaxy, a line criterion as 
described in the text on the basis of the EWs of the CaT and Mg I lines, and both criteria combined. Percentage levels are the same 
in all panels and always relative to the total number of foreground stars towards that galaxy. 



most cases either the velocity criterion is the most effective or 
performs rather comparably to the line criterion. However, for 
Canes Venatici I, Ursa Major I and Willman 1, the line criterion 
outperforms a velocity criterion. The difference is very signifi- 
cant for the latter two objects. For the three galaxies Hercules, 
Leo II and Leo IV, adding an extra criterion which is gravity- 
dependent will also pay off significantly. 

For all galaxies studied, a combination of both criteria would 
lead to an excellent cleaning of the samples, leaving only up to 
a few percent of the foreground contamination. This is partic- 
ularly important for such intrinsically faint small galaxies, that 
have a sparsely populated RGB, where the samples of targets are 
bound to be small. In such cases, even if the number of contam- 
inants is low, it can still represent a large fraction of the overall 
sample of observed stars. For example, let us take as a refer- 
ence a faint system such as UMa I, that has a l uminosity in V 
band of 1.4 ± 0.4 x 10 4 L jMarfin et al. 2008), and calculate 
what would be the expected number of RGB stars with respect 
to MW interlopers for such a system. Using the publicly avail- 
able synthetic color-magnitude diagram code provided by the 



BaSTI tearrQ and assuming that UMa I can be approximated by 
a 13 Gyr old stellar population of [Fe/H]= -2 dex, about 20-25 
stars are produced in the region of the RGB above the HB (the 
numbers are from 3 different random realizations) for a system 
of ~ 2 x 10 4 L G . For an exponential surface brightness profile, 
about 60% of the light is produced within 1 exponential radius 
- and therefore 60% of the number of stars are contained within 
1 exponential radius, assuming no mass se gregation. Using 6. 7 
arcmin as the exponential radius of UMa I (Mar tin et alj |2008). 
this would mean that approximately 12-15 RGB stars would be 
found over an area of ~140 arcmin 2 . Along the l.o.s. to UMa I, 
the Besancon model gives approximately 40 MW contaminants 
over such an area; using the velocity criterion on its own would 
result into retaining about 10 interlopers in the sample, i.e. in 
comparable number to the RGB stars of the dwarf over the same 
area. Using both a velocity and line criterion on the other hand 
would only leave 1-2 contaminants in the sample. 

Among the galaxies studied in this section, there are a few 
UFDs for which reaching about 2 mag down the HB appears to 



1 http : albione. oa-teramo inaf . i t . 
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Fig. 5. MglEW and CaT EW for the foreground stars drawn 
from the Besancon model towards the direction of the dwarf 
spheroidal galaxies Sextans (top row), Sculptor (middle row) 
and Fornax (bottom row). The left panels show all the fore- 
ground stars, while the right panels the distribution of fore- 
ground stars left after applying a velocity criterion of 3cr around 
the systemic velocity of the galaxy. The dashed line indicates 



be particularly important in order to acquire statistically signifi- 
cant number of stars, such as for UMa I, UMa II, Leo IV, Com, 
CVn II and Wil 1. We explore the performance of the method 
to fainter magnitudes in the Appendix. We restrict the analysis 
to the most compelling cases, UMa I, Leo IV, Wil 1, since for 
UMa II, Com and CVn II the velocity criterion on its own is 
already very effective. 

4. The observed trend of Mg I vs CaT EW 

We subsequently explore the observed trend of Mg I EW ver- 
sus XWcaT m actual data and how it compares to the predic- 
tions from the synthetic spectra, using the DART data-set for 
the 3 classical dSphs Sextans, Sculptor and Fornax. In this way 
we are able to empirically explore the use of the Mg I and CaT 
lines EWs for these dSphs with different properties - i.e. lo- 
cated along different line-of-sights, with different systemic ve- 
locities and whose stars cover different ranges of metallicities 
(see Table Q] for a summary). We used our intermediate resolu- 
tion spectroscopic observations of these galaxies carried out at 
the VLT/FLAMES using the GIRAFFE spectrograph in Medusa 
mode with the LR8 grating, which covers the wavelength region 
8206 -9400 A and therefore includes both the nIR CaT lines 
and the Mg I line at 8806. 8A . We refer the reader to the original 




-1.0 



Fig. 6. Comparison between MglEW measurements for stars 
with double measurements for the Sextans, Sculptor and Fornax 
dSphs (59, 137, 289 stars, respectively). The figure shows the 
distribution of Mg I EW differences from integrated flux as a 
function of S/N for the stars with S/N per A> 10 and esti- 
mated error in velocity < 5 km s _1 for each measurement. The 
weighted mean in the MglEW difference (AEWj ntj MgiX rms 
dispersion and scaled median absolute deviation from the me- 
dian are 0.01 + 0.02 A, 0.09 ± 0.03 A, 0.08 A for Sextans, 
T 0.04 + 0.03 A ,0.3 + 0.03 A , 0.19 A for Fornax, -0.02 ± 0.01 
A , 0.15 ± 0.02 A , 0.15 A for Sculptor. Assuming a similar S/N 
for the individual measurements of each star, a fit to the m.a.d. 
of these differences as a function of S/N yields "V2 x cr EW = 
V2 x 3.6/(5 IN). The solid lines indicate the 1-cr region for the 
error in the difference of EW ; the dashed lines show the relation 
used in lBattaglia et al.l(l201 ll) . The asterisks indicate the median 
AEWi ntj Ms;i per S/N bin, with the scaled m.a.d. shown as error- 
bar. 



papers for a general presentation of the data-set (fTblstov et al 
1 20041: [Battaglia et al|2006|:lHelmi et al.ll2006l:IStarkenburg et al 
bold iBaTtaglia et al l201ll) . 

The data-reduction p rocedure adopted is described in detail 
in Battaglia et al. (2008b). Here we briefly remind the reader that 
the EW of the Mg I line was derived by integrating the flux over 
6 A around the central wav elength of the line as i t was done for 
the synthetic spectra and in Batta glia et al.l (1201 ll) . For this line, 
which is smaller then the CaT lines, this EW estimator provides 
less noisy EW measurements with respect to a Gaussian fit to the 
line. 

The error in EW of the Mg I line was derived by using the 
stars with double measurements (see Fig. |6j. Assuming a simi- 
lar S/N for the individual measurements of each star, we obtain a 
typical error in the EW^ of the Mg I line, o~ ew? 0.06 A for Sextans, 
0.13 A for Fornax, 0.11 A for Sculptor as estimated from the 
scaled median absolute deviation from the median (m.a.d.) in the 
distribution of Mg IEW differences. We derive the trend of <x EW 
with S/N by fitting the value of the scaled m.a.d. of the differ- 
ence in EWnig as a function of S/N in bins containing at least 40 
stars. Since the data-sets for Sextans, Sculptor and Fornax have 
been obtained using the same instrument, observing strategy and 
data-reduction, we join the data for the 3 galaxies. We obtain 
AEW Mg = V2x<r EW = V2x3.6/(S/A0 which would correspond 
to an error on the individual measurements of <x EW = 3.6/(5' IN). 
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Figure [6] show s how this relation differs from the one used in 
iBattaglia et al.1 (1201 ll) . where only Sextans was analyzed. 

In Figure |7] we show the behavior of the Mg IEW as a func- 
tion of the line-of-sight velocity in the heliocentric system for 
the 3 dSphs. Here only spectra of higher quality (with S/N/A > 
20, velocity errors < 5 km s _1 and for which the difference be- 
tween the Mg I EW estimated from the flux integration and the 
Gaussian fit is within 5 times the scaled m.a.d. of the distribu- 
tion of Mg I |.EWint - £WGaul) are shown. This has the purpose 
of allowing us to analyze the behavior of the Mg I line without 
being affected considerably by the noise in the measurements. 
The dotted lines indicate the region of kinematic membership 
adopted in previous studies (+3cr from the systemic velocity of 
the galaxy for Sculptor and Sextans, and +2.5cr for Forna?(Q). In 
previous works stars within these velocity regions were consid- 
ered as RGB stars probable members of the dSph. 

As the figure shows, the majority of what are clearly contam- 
inant stars even by eye (i.e. those stars with line-of-sight velocity 
well outside the region of kinematic membership) can be weeded 
out on the basis of a kinematic selection. However, from the fig- 
ure it is also clear that the velocity distribution of contaminants 
extends in, and also beyond, the region of kinematic member- 
ship. Obviously, when applying a velocity criterion on its own, 
there will be contaminants with velocities consistent with mem- 
bership, and therefore classified as members. This is especially 
clear for the case of the Fornax dSph, whose low systemic veloc- 
ity falls well within the velocity distribution of MW foreground 
stars. Concerning the distribution in Mg IEW, it can be seen that 
in all 3 cases the majority of contaminants show a rather dis- 
tinct distribution in Mg I EW with respect to the bulk of proba- 
ble dSphs members, i.e. they predominantly occupy the locus of 
Mg I EW > 0.5 A . However a few differences are present among 
the galaxies, both in the distribution of probable kinematic mem- 
bers and non-members in terms of their Mg I EW. 

In the Sextans dSph 99% of the stars with velocities within 
the region of kinematic membership have Mg I EW < 0.5 A , 
which says that it is very unlikely that RGB stars - genuine 
members - will be found with MglEW > 0.5 A. Therefore an 
hard-cut on the value of the Mg I EW as applied by B 1 1 appears 
a good solution to weed out as many contaminants as possible 
with line-of-sight velocities consistent with kinematic member- 
ship, and at the same time retain in the sample the giants that are 
members of this dSph. 

The same hard-cut though cannot be applied unchanged to 
the Sculptor and Fornax dSphs. In the case of Sculptor, it is 
clear from Fig. [7] that a larger fraction of stars well outside 
the region of kinematic membership have Mg I EWs < 0.5 A 
with respect to Sextans. This is because of the different mix of 
MW stars along the line-of-sight to this galaxy (see Sect [3] and 
Figure|5j- Furthermore, about 5% of kinematic members extend 
to MglEW > 0.5 A, so that blindly imposing the same cut-off 
as Sextans could result in discarding a considerable fraction of 
genuine members of the dSph. In the case of Fornax, all but one 
of the kinematic non-members have Mg I EW > 0.5 A , but also 
the kinematic members extend to large values of the Mg I EW, 
ranging from ~0-l A. Therefore applying the same cut-off in 
Mg I EW would allow to remove all contaminants, at the price 
of selecting out about 40% of the stars with velocities consistent 
with membership! 



2 This stricter kinematic criterion was used in previous work to re- 
duce the foreground contamination given the low systemic velocity of 
Fornax 



The reason for this is clearly seen in Fig.|8]where we plot the 
distribution in the CaT YW- Mg I EW plane for highly proba- 
ble members of the dwarf and highly likely non-members. Since 
these are real data, and we do not know a priori what is a con- 
taminant and what a genuine member to the dSph, to increase the 
probability of selecting members and non-members and make 
this distinction easier using only l.o.s. velocities, we plot the 
stars within 2-cr from the systemic as members, and the stars 
beyond 4-cr from the systemic as non-members. We emphasize 
that here we are interested in the behaviour of the bulk of stars in 
both categories. For the likely members, a correlation is present 
between the 2Wc a T and the Mg I EW, so that stars with larger 
SWcaT display larger MglEW. At £W Ca T > 6 A, the typical 
Mg I EW is >0.5 A , and in Fornax, which is more metal-rich 
than Sextans and Sculptor, most of the stars in our data-set oc- 
cupy the region XWc a T ^ 6 A. The correlation in Mg I EW vs 
£WcaT observed in the data is consistent with the the synthetic 
spectra analysis of Sect. [2] and appears similar for the differ- 
ent dSphs at similar values of EW. This confirms the result of 
Sect. |2]for which a general criterion can be applied in MglEW 
vs SWcaT to separate contaminants and probable members of 
the dSphs. Obviously in the case of observations the errors in 
the measurement of the EWs of the lines need also to be taken 
into consideration when comparing the location of the star in the 
Mg I vs SWcaT plane with the line separating dwarfs and giants 
derived from the models. 

4.1. Applying the line criterion to data 

In the previous sections we have proven the applicability and 
validity of the method; it is therefore interesting to see how the 
application of the line criterion to available data-sets, e.g. DART, 
would perform. Here when referring to "contaminants" we refer 
to stars classified as such on the basis of the line criterion. 

As visual inspection of Figure [8] already suggests, a signifi- 
cant number of stars which have l.o.s. velocities consistent with 
being highly probably members of the Fornax dSphs are actu- 
ally found above the dividing line between dwarf and giant stars 
derived in Sect. [2] and therefore are most likely MW interlop- 
ers that were misclassified because went unrecognized even by 
a strict velocity selection. A smaller number of such objects is 
also present in Sextans and Sculptor. 

We apply the method to those stars that pass our usual se- 
lection criteria (i.e. S/N>10/A, velocity errors < 5 km s" 1 and 
CaT |EW int - EW g aul < 2 Afl In this work we account for the 
errors in MglEW in a simple way, i.e. we consider as "giants" 
those stars with EW Mg i lower or within 1 cr EW from the Mg I EW 
obtained by applying Eq. Q](see dividing line in Figure [3J. We 
define as "contaminants" those stars with MglEW more than 
lcr EW larger than the value obtained by applying Eq.Q] A more 
detailed approach, which would for example take into account 
the combined probability of membership from the velocity and 
Mg I EW distribution in a statistical way, is outside the scope 
of this paper; we ref er the reader to Wal ker et al.l d2009bl) and 
Mart inez et al.l d201 lh for works dealing with the issue of mem- 
bership with comprehensive statistical analyses. 

The Sextans dSphs was previously analyzed; here we just 
mention that by considering the error in Mg I line in the selec- 

3 Since part of the targets for Sculptor were chosen from preliminary 
photometry, a handful of them turned out to have colors and magnitudes 
well outside the overall selection box for RGB stars from the definitive 
photometry. We exclude these objects by restricting ourselves to V-I> 
0.5 and V > 16.5. 
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Fig. 7. MglEW versus line-of-sight velocity in the heliocentric system (squares with error-bars) for the Sextans, Sculptor and 
Fornax dSphs (from left to right). The vertical dotted lines show the region of kinematic membership adopted in previous studies 
(i.e. with velocities within ±3cr from the systemic velocity of the dSph for Sextans and Sculptor, and ±2.5cr for Fornax). 



Sextans Sculptor Fornax 




CaT EW [A] CaT EW [A] CaT EW [A] 



Fig. 8. MglEW versus SWcaT for the Sextans, Sculptor and Fornax dSphs. The filled orange squares with error-bars show the 
stars which are highly probable members (i.e. with velocities within +2cr from the systemic velocity of the dSphs) and the black 
crosses the stars which are highly likely to be MW contaminants (i.e. with velocities at least 4<x away from the systemic velocity 
of the dSphs). The solid line indicates the theoretical relation derived in Sect. [2] Note that because of the measurement errors the 
distribution of highly probable members on this plane is inflated and in some cases goes over the separation line between dwarf and 
giant stars. 



tion, only 3 stars are considered as contaminants, instead than 
the 6 stars rejected in B 1 1 . This has no influence on the results 
other than bringing back into the sample stars which had dis- 
crepant metallicities with respect to the body of the distribution, 
indication that perhaps this selection may be too generous. The 
rejected stars account for less than 2% of the overall sample. 
This may seem surprising since Sextans is a very low surface 
brightness object, with a CMD heavily contaminated by MW 
stars, and therefore one may expect a large amount of contami- 
nants. However, the systemic velocity of the object, far from the 
one of the MW disks, appears to allow for a good cleaning of the 
sample just on the basis of the l.o.s. velocity. 

Among the stars within the region of kinematic membership, 
Fornax and Sculptor both have a percentage of ~5% of stars clas- 
sified as contaminants according to the line criterion. Note that, 
for Fornax, kinematic members were chosen from a 2.5a cut in 
our previous works; if we consider the kinematic membership 
region as the velocity range between ±3<x from the systemic ve- 
locity, the fraction of interlopers increases to 6.5% of the total. 
These contaminant stars are found at all projected radii, but their 
relative contribution to the number of members becomes increas- 



ingly significant towards the outer parts, accounting for example 
for about 45% and 35% of the kinematic members at R > 1 .0 deg 
for Fornax and Sculptor, respectively. For comparison, between 
0.5 and 1.0 deg (appr ox 0.4 and 0.8 nominal tidal r adii, accord- 
ing to the value from llrwin & Hatz idimitrioul l 19951) the fraction 
drops to 6.3% in Sculptor. The fraction of contaminants found 
for these two objects are very similar to each other, even though 
Table 1 would have suggested that a much larger fraction would 
be found for Fornax. This is most likely due to a combination 
of factors such as the different area covered and the much larger 
surface brightness of Fornax with respect of Sculptor, which per- 
mits to pick out a larger number of member stars from the fore- 
ground. 

Figure [9] (top) shows the distribution of the analyzed stars 
for Sculptor and Fornax on the CaT [Fe/H] versus heliocentric 
l.o.s. velocity plane. The combination of these quantities can al- 
ready provide information on the probability of membership of a 
star to the MW or to the dwarf galaxy, as on this plane the great 
majority of those stars which are clearly unbound to the dwarf 
galaxy (MW stars) has a different distribution than the velocity 
members (likely RGB stars of the dSphs). For example, it is clear 
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Fig. 9. Application of the line method to the Sculptor (left) and Fornax (right) dSphs, using the VLT/FLAMES DART data-set. 
Top panels: Distribution on the CaT [Fe/H] versus l.o.s. heliocentric velocity for the stars with S/N/A >10, error in velocity < 5 
km s _1 and CaT |EW; nt - EW gau | < 2 A. The black squares with error bars show those stars whose Mg IEW was either lower or 
within lcr EW from the Mg IEW given by Eq. Q] while cyan (grey in the black and white version) squares with error bars show 
those stars whose Mg IEW was more than l<x EW larger than the Mg IEW given by Eq. [1] The vertical lines indicate the region of 
kinematic membership from our previous works. Bottom panel: CaT [Fe/H] versus elliptical radius for kinematic members (black: 
with Mg I EW within 1 cr EW from Eq. [T] cy an: with Mg I EW more than 1 cr EW away from Eq. Q] 



that the region with -2 < [Fe/H]< -1 is expected to be the most 
contaminated one. The figure also shows that most of the stars 
excluded by the line criterion alone indeed have CaT [Fe/H] val- 
ues in the range where most contaminants were expected, high- 
lighting the fact that the method here presented, which makes 
no a priori selection in metallicity and velocity, is indeed able to 
find interlopers. 

We note that the CaT [Fe/H] relation is empirically calibrated 
for RGB stars (above the HB) and relies on the assumption that 
all the stars can be considered at the same distance. While this 
is a reasonable assumption for the stars belonging to the dSphs 
here analyzed, this is not the case for the bulk of MW interlopers, 
which are located over a range of distances - mostly within a few 
kpc from the Sun - and have markedly different gravities than 
RGB stars. Therefore, the CaT [Fe/H] values for the stars which 
are probable non-members are not necessarily indicative of their 
true metallicity. 

Figure [9] (top) also illustrates what would be the result of 
relaxing the kinematic criterion to Act for the 2 galaxiefl: the 
fraction of interlopers becomes about 8% for Fornax and 6.5% 
for Sculptor; and the new contaminants are found at velocities 



4 Here we use the velocity dispersion values from Table 1 . 



lower than the systemic, towards the velocities of the disk, as 
it would be expected. The bottom panels of Figure [9] show the 
metallicity distributions as a function of projected radius result- 
ing from applying the line criterion in addition to previous kine- 
matic criteria: the conclusions from our previous works in terms 
of changes of metallicity with radius holds. Also the smaller ve- 
locity dispersion of metal-rich stars with respect to metal-poor 
stars is still present. For Sculptor it can be seen that several of 
the stars with [Fe/H]> -1.7 dex at R>0.6 deg may actually not 
belong to Sculptor, making the variation of metallicity properties 
with projected radius even clearer. 



One can therefore conclude that the main properties of sys- 
tems for which samples of several hundreds individual stars are 
available, such as the majority of classical dSphs, are robust. 
It is also clear that methods of the kind here proposed can po- 
tentially be very helpful for those analyses whose conclusions 
heavily rely on the properties of stars in the outer parts of these 
systems, where the ratio of unrecognized interlopers to member 
stars increases, as well for studies of several of the fainter dwarf 
galaxies such as the UFDs. 
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5. Discussion and conclusions 

In this work we have addressed the problem of identifying MW 
stars in samples of kinematically selected members stars in 
Local Group galaxies. The presence of such interlopers may con- 
taminate the derived properties of the target galaxy and this is a 
general issue in studies of resolved stellar populations of Local 
Group galaxies. 

The method that we propose is mainly meant for spectro- 
scopic samples in the CaT region whose targets are selected from 
broad band photometry to lie on the a region of the CMD cov- 
ering the RGB of the target galaxy (above or approximately one 
magnitude below the HB). This is the case for the great major- 
ity of data for "classical" dSphs and more distant Local Group 
dwarf galaxies; for several UFDs a significant fraction of the tar- 
gets lie in that magnitude range. We have explored the combined 
use of the Mg I EW at 8806.8 A and the nIR CaT 2W as a way 
of distinguishing between contaminants and members on the ba- 
sis of physical characteristics of the star, mainly gravity, since 
the targeted stars are RGB stars while the MW interlopers will 
mostly consist of dwarf stars. 

For this we have used synthetic spectra over a range of metal- 
licity, temperature and gravity covering the expected age and 
metallicity range of stars in Local Group dwarf galaxies. We 
found that a relation between the Mg I EW and the CaT ZW can 
be applied to distinguish between dwarf and giant stars above 
the HB for metallicities [Fe/H] > -2 dex and for (sub)giant stars 
below the HB for -2< [Fe/H] <- 1 . Since at lower metallicities the 
main contribution of contaminants comes from the MW stellar 
halo, which is the least dominant component in terms of stellar 
mass (and therefore number of stars), only a negligible amount 
of interlopers is anyway expected below this metallicity. 

The relation between the MglEW and the CaT 2 W de- 
rived from the analysis of synthetic spectra applies well to 
the trend observed in actual data for a sample of 3 classical 
dSphs (Sextans, Fornax and Sculptor), whose stars cover differ- 
ent ranges of age and metallicity. The line-of-sight to the target 
object and the velocity of the center of mass of the galaxy (the 
systemic velocity) are important parameters because they deter- 
mine the amount and mix of contaminants that will be present in 
the kinematically selected samples of stars. This means that there 
will be objects for which the l.o.s. velocity selection is already 
very efficient, while other ones that greatly benefit from an ad- 
ditional criterion such as the one here explored. We find that the 
addition of this method to the information from the l.o.s. veloc- 
ity can be particularly beneficial for example for Canes Venatici 
I, Leo II, Fornax, Hercules, Ursa Major I and Willman 1. For the 
latter four galaxies, the use of the Mg I EW method on its own is 
considerably more efficient than the selection in l.o.s. velocity. 

It should be noted that the contamination due to giant stars 
from the MW is not removed by this method, the proposed rela- 
tion will retain as many giant stars as possible and throw away 
only dwarf contaminants. Another very important parameter in 
determining the ratio between contaminants and stars member 
to the dwarf is the ratio between the projected number density of 
MW and dwarf galaxy stars. Besides varying from object to ob- 
ject, this ratio changes within the object itself simply because of 
the declining surface density of stars in dwarf galaxies, so that 
weeding out interlopers will become particularly important in 
the outer parts of the systems, and for intrinsically faint objects. 

This method has the advantage of being very easily imple- 
mented, as it relies on measuring EWs for relatively large lines 
such as the nIR CaT and Mg I . These lines are close in wave- 



length, so that several existing studies of MW dSphs covering 
this wavelength region could benefit from this extra information. 
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Appendix 

Here we explore how the line criterion would perform down to 2 
mag below the HB for Uma I, Wil 1 and Leo IV, i.e. those UFDs 
for which the application of the line criterion would considerably 
improve the elimination of foreground stars but whose faintness 
makes it very difficult to acquire large number of targets even 
reaching down 1 mag below the HB. 

We do so by increasing the temperature range and includ- 
ing models as hot as 5500 K; while this results into colors only 
milder bluer than previously considered, i.e. V-I ~ 0.8, it allows 
to extend the magnitude range of the sub-giant branch boxes of 
about 1 mag at the fainter limit. 

Fig.flOlshows that qualitatively the results of Sect.[2]continue 
holding: at [Fe/H]=-2 and -1 it is possible to make a distinc- 
tion between dwarfs versus giants/sub-giant stars; at [Fe/H]= 
sub-giants overlap with both the locus of dwarfs and giants, but 
dwarfs and giants remain distinct; as before, at [Fe/H]=-3 the 
various types of stars overlap in the SWc a T - EWjyig plane. 

Obviously the actual efficiency in removing contaminants 
will change from what considered in Sect. 13.21 as by increas- 
ing the temperature range at the higher temperature end, we are 
not only increasing the number of contaminants but also includ- 
ing hotter MW stars. Table [2] shows the results when selecting 
stars from the Besancon model to have 0.8 <V-I< 2.1 and mag- 
nitude from the tip of the RGB (as done in Sect. [3} down to 2 
mag below the HB: the fraction of contaminant stars retained 
in the sample by both velocity and line criterion increases, but 
the combined efficiency is still very good, with only 3-5% of in- 
terlopers retained. It should be also noted that for Ursa Major I 
and Willman 1 the line criterion still outperforms the velocity 
criterion. 
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Fig. 10. As in Fig.[3]but imposing an upper limit of 5500 K for the selection in temperature, rather than 5250 K. 

Table 2. Properties of the foreground contamination towards the MW satellites Leo IV, UMa I and Wil 1 . The columns list: the name 
of the galaxy (1), the number of expected MW contaminants within a solid angle of 4 deg 2 along the line-of-sight to the galaxy 
(2); the percentage of contaminants retained in the sample when using only the velocity criterion (3), the line criterion (4), when 
applying both (5). 



Galaxy 


N co „, 
/ 4 deg 2 


%(N„„„) 
vel. crit. (3cr) 


%(N«„„) 
line crit. 


%(N«*) 
both (3cr) 


Leo IV 


8063 


13% 


25% 


5% 


Ursa Major I 


7147 


34% 


17% 


5% 


Willman 1 


3806 


36% 


15% 


3% 
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